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Abstract 

For direct dark matter detections, to extract useful information about the fundamental interaction 
from data, it is crucial to properly determine the nuclear form factor. The form factor for spin- 
independent cross section of collisions between dark matter particle and nucleus is thoroughly 
studied by many authors. When the analysis was carried out, the nuclei are always supposed to 
be spherically symmetric. In this work, we investigate the effects of deformation of nuclei from a 
spherical shape to an elliptical shape on the form factor. Our results indicate that as long as the 
ellipticity is not too large, such effects cannot cause any substantial effects, especially as the nuclei 
are randomly orientated in a room temperature circumstance one can completely neglect them. 



PACS numbers: 21.60-n, 24.10. Cn, 95.35+d 
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I. INTRODUCTION 



With serious astronomical observation of several decades, existence of dark matter is 
no longer doubtful. On another aspect, we definitely know that in the zoo of the standard 
model (SM) we do not have any candidates for dark matter (DM). Question is what the 
Dark Matter particles are. There are many models proposed in literature l|-5|, but unless 
they are caught by our detectors in the terrestrial laboratories or satellites 6|-l8[, one still 
cannot surely identify them. Much efforts have been made to discover the dark matter flux 
from outer space. 

Comparing with the spin- dependent cross section, the spin-independent cross section 
of dark matter particle with nucleus is much larger due to the enhancement where A is 
the atomic mass number of the nucleus as the detection material [9-15|. Even so, the cross 
sections of elastic scattering between DM particles and nuclei are still small, the present ex- 
periments have already reached 10~^^ cm^. Because of the advantage of the spin-independent 
scattering whose cross sections are larger and the theoretical treatment is relatively simpler 
than that for spin-dependent processes, nowadays, the priority of research is given to the 
study on spin-independent elastic DM-nucleus reactions. 

Since the kinetic energy of the DM particle is rather low at order of a few tens of keV, the 
impact of DM particle on nucleus is almost impossible to cause inelastic processes, thus all 
observational signals are related to the recoil of nucleus after the collision. Namely, even the 
collision occurs between DM and quarks (seldom gluons in the case), all the absorbed energy 
is passed to the nucleus to make it to move as a whole object, the recoil which may induce 
thermal, electric and light signals. For the spin-independent cross section, the particle- 
physics and nuclear-physics contributions can be separated, namely the nuclear effects can 
be factored out and included in a form factor F{q). For spherically symmetric nuclei, F{q) 
only depends on |q|. The spherical symmetry means that the nuclei are of full-shell or close 
to full-shell structures, but for most of the nuclei which are taken as the detection materials 
the shells are not completely filled out. A careful study on form factors for the not-full-shell 
structure nuclei would be helpful for extracting information about fundamental interactions 
from the data. For that case, -F(q) is not only a function of |q|, but also cos^ while the 
azimuthal symmetry is assumed. We will write it as F[q, cos6) where q = |q|. 

Nucleus is a complex many-body system, therefore extraction from data requires a 
thorough analysis on the nuclear structure. The form factors for spherical nuclei have been 
carefully studied by many authors and the results can be applied to analysis of data. In 
this work, we are going to investigate the effects of deformation of nuclei on the form factor, 
namely we will derive the form factors corresponding to the deformed nuclei with relatively 
smaller ellipticity. 

We employ several models to calculate the form factors F{q, cos 6) for nuclei with small 
ellipticity. We will take Xe and Ge which are commonly adopted as the detection materials 
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as examples to illustrate the effects of deformation. 

The paper is organized as follows, after this introduction, we present the expressions 
of the form factors derived from different models for nuclear density, and then present our 
numerical results via several figures. The last section is devoted to our conclusion and 
discussion. 



II. THE FORM FACTOR RELATED TO DEFORMED NUCLEUS 



Obviously it is reasonable to assume that the nucleus with a larger A may be only 
quadruply deformed, namely it is deformed from spherical form to ellipsoidal. In the spherical 
coordinates, the nuclear density of a nucleus with an elliptical form should be p{r, 9) which 
is a function of both radius r and polar angle 9 and the corresponding form factor should be 
written as follows, it needs to be noted that we would set ipi = ip2 in practical calculation 
for simplifying the integration. 



F{q,9,) = ^ j p{r,9,)e^'^-^^d\ 

-| PIT p27V poo 

M Jo Jo Jo 

Even though this work is aiming to find the effects of deformation of nuclei on the form 
factor, the deviation from the spherical form for the nuclei under investigation is not severe, 
therefore, we can always start from a spherical form and then make reasonable modifications 
or extension. 



A. Extension of the Two-Parameter Fermi Distribution(E2PF) 



A number of models have been proposed 16|, ll7|] to describe the nuclear charge density 



or mass density. Among them the two-Parameter Fermi distribution(2PF) is one of the 
simple models. For a spherical form the density is written as 



Po 

1 + exp( 



PW = . , .'Lr-c, , (2) 



where po is equal to 2p(r) at r = c, and z is the diffusivity of the surface. It would be 
convenient for later use to derive the mean square root radius R for the spherical 2PF model 

as 

p2PF l ^^irr'P^ /3 7 
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For an elliptical nucleus with an axial symmetry, the nuclear density in the two 



parameter Fermi distribution (E2PF) model should be extended as 18l425l|: 

pA_ 

1 + exp ^ 



Pir,0) = ^^r-^, (4) 



where 

ci9) = co{l + f32Y2o{9)). (5) 

The parameter which corresponds to the ellipticity of the nucleus characterizing its de- 
formation from a spherical form, is a small quantity for the nucleus which we concern. For 



the priori assumption of small deformation, we only keep the multiple terms up to Y20 (25|]. 

The parameter cq = I.IA^/^ and pg can be obtained by the normalization condition 
i.e. requiring the integration over whole coordinate space to be equal to the nuclear mass 
number (or total charge Ze) which is priori set for various nuclei. P2 can be obtained from 



the data book 26[, and it is —0.113 and —0.224 for ^^^Xe and ^^Ge respectively, z denotes 



the surface diffuseness. Here we choose the normalization as follows: 

p(r, e)d\ = M. (6) 




radius(fm) 



FIG. 1: Nuclear density of Xe for the Extended 2PF model. It shows the change of density with 
increasing angle from 10° to 90°, the Short Dash Dot line(green) corresponds to the case of the 
spherical 2PF model 
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momentum transfer q(1/fm) 



FIG. 2: ^^^Xe form factors for deformed nucleus Extended 2PF(E2PF) model from different direc- 
tions: 10", 30°, 45°, 60°, 90° 



In Fig.(IT]) we show the density distribution for ^^^Xe in the E2PF model. It is observed 
that from the center of the nucleus to about three fermis, the density remains unchanged in 
all directions. Then the angular distribution of the density begins to be apart for different 
angles beyond three fermis. The short dash dot (green) line is the 2PF density model when 
the nucleus is assumed to be spherical. Fig.(|2]) is the corresponding form factors, which 
are calculated by taking a Fourier transformation to the deformed nuclear density in the 
configuration space. 



B. Extension of the Folding model(EF) 



There is another commonly adopted model which is rather simple, i.e. the nucleons are 
postulated to be uniformly distributed in a sphere with a certain boundary radius. For an 
axially symmetric ellipsoidal shape, one should extend the density for a spherical form. The 
surface equation of an ellipsoid is 

- + ^ + - = 1 (7) 
or in the spherical coordinate system, it is written as 
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In an approximation, if we only keep the multiple terms to quadrupole, we can re-parametrize 
the surface equation to a more convenient one 

R{e)=Ro{l+/32Y2o). (9) 

Extending the Folding model, we set the nuclear density to be uniform inside the ellip- 
soid with radius R{0) 

Po(r,^) = ^e(r-i?(e)), (10) 



where G is the step function. Following the literature 27|], we introduce a smearing function 
pi to take care of the soft edge effect of the nucleus: 



^2 



PiW = 77^33^ exp(^), (11) 



then one should convolve po and pi to get the nuclear density 



p{r,6) = J po{r')pi{f — r')(fr' 

= j Po{r')pi{f -r')r'^dr' sin 9' d9'd(f' 

d^' sinW / J±i^e(r' - i?(e')) X 
Jo Jo 47ra^6 



(27rs2)3/2 

(r^ + r'^ — 2rr' (sin ^ sin ^' COS (</) — + COS ^' COS 6')) i^it 
exp( ^-^ )r dr . (12) 



The semi- axes a and h are set as 



a = R{d = ^)= i?o(l + /32F2o(^)) = Ro{l - \ T^/32) 
2 2 V lovr 



h = i?(^ = 0) = /?o(l + /32l^2o(0)) =/2o(1 + 2a/-|-/32). (13) 

V Idtt 

In the extended Folding (EF) model, we can also calculate the mean square root radius 
R^^ . For a spherical nucleus, one may equate the mean square root radius obtained in the 
2PF and Fold models, thus he acquires the spherical radius Rq for the Fold model 



16 



28 



Ro = ^c^ + ^7r2a2 - 5s2, c ~ (1.23A^/^ - 0.6)fm, s = 0.9fm, a = 0.52fm. (14) 

As aforementioned, the deformation makes the shape of the nucleus slightly deviate 
from a spherical form, we can still use the above relation achieved for spherical nuclei and 
set Rq to be the parameter in Eq.([9]). 
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A Fourier transformation would bring the nuclear density to the expected form factor 
6). It is noted that now the form factor is also direction-dependent. 

F{q) = y"po(r')pi(r-r')rfVV'i"rfV 

= j po{r')d\' j pi(r-r')ciW'i"ciV 

noo POO 

o(r')t^V / pi{u)e''^''d\ (setting r - r' = u) 
Jo 



Pol 

Jo 

/*00 POO 

Jo Jo 

POO POD 



= Fo(g)Fi(g), (15) 

Fi{q) = j pi(r)e''i"rfV 

= e-^"/'. (16) 

We take a trick to make the integration easier as r is described in the cylindrical coordinate 
while q is described in the spherical coordinate as: 

X = t cos (pi \ \ Qx = Q sin 62 cos (p2 
y = tsin(pi > <^ Qy = q sin 62 sin 



z = z 



qz = q cos 62 



then 

q ■ r = t cos (piq sin 6'2 cos 02 + t sin 0ig sin 62 sin </)2 + cos 6*2, 

where t = \J + y"^ = y/r'^ — z'^. 
Thus we obtain 

Foiq,02) = J po{r,e)e'^-^d'r 

Ana^b 
Ana^b J 

The parameters a and 6 are the semi-axes defined above. Thus the form factor in the EF 
model can be written as: 

Fiq,e2) = Fo{q,e2)F,iq). (18) 
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In Figure [3l the ^^^Xe density distribution determined by the EF model is shown, while 
the corresponding form factors are given in Fig. 4. The short dash dot (green) line corresponds 
to the spherical form of the nucleus. For the spherical nucleus, it has already been known 
that the form factor obtained with the 2PF model is very close to that determined by the 
Folding model[l3]. Thus we will also make a comparison between the E2PF and EF form 
factors at the end of the paper. Then we will present the third model to do the same job in 
the following section. 



0.17- deformd nuclear density of Extended Folding model (EF) 




radius(fm) 



FIG. 3: ^^^Xe form factors for deformed nucleus of Extended 2PF(E2PF)model from different 
directions: 10^ 30°, 45°, 60°, 90° 



C. The Nilsson Mean Field(NMF) 

Above, we use two simplified models (E2PF and EF) to derive the form factors for 
deformed nuclei. The advantage is that the models are simple and we can obtain analytical 
solution which is convenient for illustrating the characteristics of the form factors, but might 
be too simplified. Now we turn to use a more realistic model. 

In this subsection, the form factors for deformed nuclei are obtained in the Nilsson 
modified oscillator model, then using ^^^Xe as an example, we present the results in some 
figures. 

Below, let us briefly review the model and show how we apply it to study the concerned 
form factor. 
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' Xe form factors for deformed Extended Folding(EF) 
model from different angle 




— -90° 

— 60° 

45° 

■■■■30° 

10° 

-- spherical 

spherical nucleus form factor of Fold" 
model(green line) 



0.2 0.3 0.4 0.5 0.6 0.7 O.f 

momentum transfer q(1/fm) 



0.9 1.0 



FIG. 4: ^"^^Xe form factors for deformed nucleus determined in the EF model for different directions: 
10°, 30°, 45°, 60°, 90° 



In the Hamiltonian of the Nilsson model the potential for an axially symmetric harmonic 



oscillator can be written as 



H = + -M[ul{x' + y') + uy] - Cs ■ 1 - Dl', 



(19) 



where Cs ■ 1 is the spin-orbit coupling, and Dl^ flattens the bottom of the potential. 

A deformation parameter 6 is introduced to reflect the axial symmetry for the deformed 
nuclei as 



The equipotential surface encloses a constant volume if 



uj^ujyijJz = const. 



Then we have 



Uo[l 



16 
— ( 
27 



31 1/6 



1^00- 



(20) 
(21) 

(22) 
(23) 
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The Hamiltonian can be decomposed into three pieces as 

H = Hsp + H,^-Cs-l- D\\ (24) 

Hsp = — + -mu^r^, (25) 

H,, = -muy\e2P2{cosd). (26) 

It would be convenient to use dimensionless coordinates and parameters which are 
defined as 



P = ^^r, (27) 

C = 2iihujoo, (28) 

D = 2^/^ ~ nhjoQfi. (29) 
Then the Nilsson Hamihonian can be further written as 

H = hwoiHo - ^esPs) - ^/iWoops ■ 1 + fi{f- < >n)1 (30) 

Ho = l{-Vl + p% (31) 

where < P >= A^(A^ + 3)/2 is an average over all states within the iV— th shell, and hu^Q ~ 
41A-1/3 MeV. 

If the octupole and hexadecupole deformations are considered, the Hamiltonian would 
become more complicated as 

H = huo{Ho + p2(_2/3e2P2 + esPs + e^P^)) - KhuJoo[2l ■ s + fi{l^- < 1^ >;v)] (32) 

and it is the Hamiltonian we are going to use in the later part of this paper. 

The Nilsson wavefunction is constructed with the spherical harmonic oscillator basis 

\Nijn >, 

= ^a;,cl|0>, (33) 

a 

where Ua is a coefficient, a refers to a set of quantum numbers (njlfl) of the harmonic- 
oscillators. 

The nuclear density thus is expressed as 

p = 5^(vl/t(7r)vl/,(7r) + vl/t(vr)vl;,(7r)) + ^(vl>t(^)vl;,(z.) + vl>t(^)vl>,(^)), (34) 

1=1 i=l 

where i represents the time-reversed states. 

In this paper, the major shells under consideration are from to 9 for proton and 
neutron, respectively. The quadrupole, octupole and hexadecapole deformation parameters 
are determined by experiments 32f], which are -0.108, and 0.027, respectively. 
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momentum transfer q(1/fm) radius(fm} 



FIG. 5: The right panel shows the dependence of the ^ Xe density on the directions from 10° to 
90° obtained in the Nilsson Mean Field model, the left panel is the form factor 



Performing a Fourier transformation on the nuclear density, we obtain the concerned 
form factor. On the right panel of Fig. ([5]), the angular dependence of the nuclear density 
of ^^^Xe is shown. Fig. ([5]) plots the NMF form factors F{q,6) with various angles. The left 
panel of Fig.(j6]) compares the form factors at direction of ^ = 7r/6 obtained with the three 
models: the EF, E2PF, and NMF, whereas the right panel is the corresponding densities. 
Fig-© shows difference of the form factors F(g, 7r/6) for ^'^Ge and ^'^^Xe, as well as their 
density distributions. 




FIG. 6: The form factors F(q,7r/6) and density p{r,ir/6) obtained in three different models: 
E2PF,EF,NMF for ^^iXe 
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momentum transfer(1/fm) radius(fm} 

FIG. 7: the left graph is the form factors F(q,7r/6) for the '^'^Ge and ^^-"^Xe with NMF model, the 
right graph is the density distributions p(r, vr/G). 

III. SUMMARY 

The aim of this work is to study if a small deformation of nuclei can induce observable 
effects on the form factors for the direct detection of dark matter flux. The form factor, 
no matter the nuclei are of spherical or deformed shapes, say ellipsoidal, must satisfy two 
normalization conditions. First, the nuclear density must be normalized as 

j p{r, 9, (p)d^r = M, (35) 

where M is the total mass of the nucleus. It is independent of the shape of the nucleus. 
Secondly, the form factor must satisfy the condition 

F(|q| = 0) = F(0) = l. (36) 

This condition does not depend on polar and azimuthal angles. With these two conditions, 
one can adopt different models for the nuclear density and then carry out a Fourier trans- 
formation to convert the nuclear density from the configuration space into the momentum 
space to gain the form factor which corresponds to non-zero momentum transfer. 

In this work, we start with the spherical nuclei and adopt three models which are 
commonly employed to study the nuclear effects. Then we extend them to deformed shapes 
by including polar angle dependence in the density while an axial symmetry is assumed for 
simplicity. With the three models we obtain the form factor for the nuclei whose shape 
slightly deviates from spherical form, namely their ellipticity is relatively small. 

We notice from the figures shown in the text that the form factors are not very much 
apart from that for spherical form, indeed the dependence of the form factor on |q| for 
6' = 7r/4 is rather close to that for spherical shape. 
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Especially, if there is not a strong magnetic field to polarize the nuclei at very low 
temperature, the nuclei in the detection material would be randomly oriented, the polar 
and azimuthal angles would be averaged and the deformation effects should be eventually 
smeared out. 

Therefore, our conclusion is that unless one can keep the detector at very low tem- 
perature such as the CDMS detector and apply a strong magnetic field to it, the effects of 
deformation of nuclei can be safely ignored. 
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